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Abstract. Say that an edge of a graph G dominates itself and every 
other edge adjacent to it. An edge dominating set of a graph G = (V, E) 
is a subset of edges E' C E which dominates all edges of G. In particular, 
if every edge of G is dominated by exactly one edge of E' then E' is a 
dominating induced matching. It is known that not every graph admits 
a dominating induced matching, while the problem to decide if it does 
admit it is NP-complete. In this paper we consider the problems of finding 
a minimum weighted dominating induced matching, if any, and counting 
the number of dominating induced matchings of a graph with weighted 
edges. We describe an exact algorithm for general graphs that runs in 
O* (1.1939™) time and polynomial (linear) space. This improves over any 
existing exact algorithm for the problems in consideration. 



Keywords: exact algorithms, dominating induced matchings, branch 
& reduce 



1 Introduction 



Under the widely accepted assumption that P ^ NP there are several 
problems with important applications for which no polynomial algorithm 
exists. The need to get an exact solution for many of those problems 
has lead to a growing interest in the area of design and analysis of exact 
exponential time algorithms for NP-Hard problems j!3!24j . Even a slight 
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improvement of the base of the exponential running time may increase 
the size of the instances being tractable. There has been many new and 
promising advances in recent years towards this direction |1|2| . 

In this paper we give an exact algorithm for the weighted and count- 
ing version of the NP-Hard problem Dominating Induced Matching (also 
known as DIM or Efficient Edge Domination) which has been extensively 
studied [9 8 16 19 2 0|4|5|3|7j . Further notes about this problem and some 
applications related to encoding theory, network routing and resource al- 
location can be found in |14|18j . 

The unweighted version of the dominating induced matching problem 
is known to be NP-complete |14j . even for planar bipartite graphs of 
maximum degree 3 [3] or regular graphs [9]. There are polynomial time 
algorithms for some classes, such as chordal graphs |19j . generalized series- 
parallel graphs [19] (both for the weighted problem), claw-free graphs [7], 
graphs with bounded clique- width [7], hole-free graphs [3J, convex graphs 
[IB] , dually-chordal graphs [4j, JV-free graphs [5], bipartite permutation 
graphs [20] (see also [6]). 

If P 7^ NP it is not possible to solve this problem in polynomial time, 
hence it becomes important to improve the exponential algorithm in order 
to identify the instances that can be solved within reasonable time. 

A straightforward brute-force algorithm to solve weighted DIM can 
be achieved in 0*(2 n ) time and polynomial space. 

The minimum weighted DIM problem can be expressed as an instance 
of the maximum weighted independent set problem on the square of the 
line graph L(G) of G, and also as an instance of the minimum weighted 
dominating set problem on L(G), by slightly way described in |4|21j for 
unweighted DIM problem. 

The minimum weighted dominating set can be solved in O*(1.5780 ra ) 
time |12j . while the maximum weight independent set can be solved in 
0*(1.4423 n ) time by enumeration of all maximal independent sets |23j. 
To the best of our knowledge there are no better method to obtain the 
maximum weighted independent set (a better algorithm O*(1.2209 n ) for 
unweighted maximum independent set is due by |llj). Hence the DIM 
problem for a graph G can be solved by using this algorithm in L 2 (G), 
which runs in 0*(1.4423 m ) time. 

For the minimum weighted DIM this algorithm behaves better than 
the brute-force alternative whenever 1.4423" 1 < 2 n , this is, m < 1.8926n. 

The paper [17] shows how to solve the DIM problem in 0*(1.7818 n ) 
time and polynomial space while the same algorithm runs in 0(n + m) 
time if the graph has a fixed dominating set. In the same work another 
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approach based on enumerating maximal independent sets was developed 
and allows to solve both DIM problems (minimum weighted problem 
and counting problem) in 0*(1.4423 n ) time and polynomial space. For 
the counting problem, there exists algorithms such as [10] which can be 
used to count the number of MWIS's in 0*(1.3247 n ) time, leading an 
0*(1.3247 m ) time and polynomial space algorithm to count the numbers 
of DIM's. 

Comparing with the straightforward brute-force algorithm, it is con- 
venient to use it as long as 1.3247 m < 2 ra , and this occurs whenever 
m < 2.4650n. 

There are NP-complete instances of the DIM problem where the num- 
ber of edges in G is relatively low such as for planar bipartite graphs of 
maximum degree 3 [3], where m < 1.5n. Therefore using transformations 
and exact algorithms for MWIS or for counting MWIS's is better than 
using the brute-force algorithm. Note however that cases where brute- 
force algorithm is not convenient strongly relies on the number edges. 
For instance, for a graph with 0(n) edges such that m ~ 3n the MWIS 
algorithm behaves better than the brute-force one. 

In this paper, we propose an algorithm for solving the weighted DIM 
problem and the counting DIM's problem, in 0(m-1.1939 n ) € 0*(1.1939 n ) 
time and 0(m) space in general graphs which improves over the existing 
algorithms for these problems. We employ techniques described in |13| for 
the analysis of our algorithm, and as such we use their terminology. 

The proposed algorithm was designed using the branching & reduce 
paradigm. More information about this design technique as well as the 
running time analysis for this kind of algorithms can be found in jT3]. 

2 Preliminaries 

By G(V,E) we denote a simple undirected graph with vertex set V and 
edge set E, n = |V| and m = \E\. We consider G as a weighted graph, that 
is, one in which there is a non-negative real value, denoted weight(vw) 
assigned to each edge vw of G. If v £ V and V' C V, then denote by 
N(v), the set of vertices adjacent (neighbors) to v, denote d(v) = \N(v)\ 
the degree of the vertex, denote by G[V] the subgraph of G induced by 
V , and write Nyi(v) = N(v) D V . Some special graphs or vertices are 
of interest for our purposes. A graph formed by two triangles having a 
common edge is called a diamond. By removing an edge incident to a 
vertex of degree 2 of a diamond, we obtain a paw. Finally, a vertex of 
degree 1 is called pendant. 
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Given an edge e £ E, say that e dominates itself and every edge 
sharing a vertex with e. Subset E' C £7 is an induced matching of G if each 
edge of G is dominated by at most one edge in E'. A dominating induced 
matching (DIM) of G is a subset of edges which is both dominating and 
an induced matching. Not every graph admits a DIM, and the problem of 
determining whether a graph admits it is also known in the literature as 
efficient edge domination problem. The weighted version of DIM problem 
is to find a DIM such that the sum of weights of its edges is minimum 
among all DIM's, if any. The counting version of the problem consists 
on counting the amount of DIM's a graph has. It is easy to see that the 
weighted version and the counting version of the problem are harder than 
the unweighted one. If the graph G has negative weights the problem 
can be solved using the same algorithm that solves the problem for non- 
negative weights. Let —M be the minimum weight among all edges of G, 
modify the weights of G by adding M to the weights of all edges. It is not 
hard to see that every DIM is a maximum induced matching, and hence 
the number of edges of every DIM in G is the same. Therefore the optimal 
solution for the modified graph is the same that the optimal solution for 
the original graph. 

We assume the graph G to be connected, otherwise, the DIM of G is 
the union of the DIM's of its connected component, and so we can restrict 
to the connected case. 

We will use an alternative definition [Sj of the problem of finding a 
dominating induced matching. It asks to determine if the vertex set of 
a graph G admits a partition into two subsets. The vertices of the first 
subset are called white and induce an independent set of the graph, while 
those of the second subset are named black and induce an 1-regular graph. 

A straightforward brute-force algorithm for finding the DIM of a graph 
G consists in finding all bipartitions of V(G), color one of the parts as 
white, the other as black, and checking if the result is a valid DIM. The 
complexity of this algorithm is 0(2 n ■ m) € 0*(2 n ). 

3 Extensions of Colorings 

Assigning one of the two possible colors, white or black, to vertices of G 
is called a coloring of G. A coloring is partial if only part of the vertices 
of G have been assigned colors, otherwise it is total. A black vertex is 
single if it has no black neighbor, and is paired if it has exactly one black 
neighbor. Each coloring, partial or total, can be valid or invalid. 
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Next, we describe the natural conditions for determining if a coloring 
is valid or invalid. 



Definition 1. ; RULES FOR VALIDATING COLORINGS: 

The following are necessary and sufficient conditions for a coloring to be 

valid: 

A partial coloring is valid whenever: 
VI. No two white vertices are adjacent, and 

V2. Each black vertex is either single or paired. Each single vertex has 
some uncolored neighbor. 

A total coloring is valid whenever: 

V3, No two white vertices are adjacent, and 
V4. Each black vertex is paired. 

Lemma 1. There is a one-to-one correspondence between total valid col- 
orings and dominating induced matchings of a graph. 

Proof: It follows from the definitions, a 

For a coloring C of the vertices of G, denote by C -1 (white) and 
C -1 (black), the subsets of vertices colored white and black. A coloring 
C' is an extension of a C if C^ 1 (black) C C'^ 1 (black) and C' 1 (white) C 
C'- l (white). For V , V" C V(G) if C is obtained from C by adding to 
it the vertices of V' with the color black and those of V" with the color 
white then write C = C' U BLACK (V) U WHITE(V"). Note that a 
total valid coloring can be only an extension of partial valid colorings and 
itself. 

Given a partial coloring C, the basic idea of the algorithm is to iter- 
atively find extensions C of C, until eventually a total valid coloring is 
reached. It follows from the validation rules that if C is invalid, so is C . 
Therefore, the algorithm keeps checking for validation, and would discard 
an extension whenever it becomes invalid. 

Basically, there are two different ways of possibly extending a coloring. 
First, there are partial colorings C which force the colors of some of the 
so far uncolored vertices, leading to an extension C of C. In this case, 
say that C has been obtained from C by propagation. The following is 
a convenient set of rules, whose application may extend C, in the above 
described way. 
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Lemma 2. ; RULES FOR PROPAGATING COLORS: 



The following are forced colorings for the extensions of a partial col- 
oring of G. 

PI. In an induced diamond, degree-3 vertices must be black and the re- 
maining ones must be white 
P2. The neighbor of a pendant vertex must be black 
P3. Each neighbor of a white vertex must be black 

P4. Except for its pair, the neighbors of a paired (black) vertex must be 
white 

P5. Each vertex with two black neighbors must be white 

P6. If a single black vertex has exactly one uncolored neighbor then this 

neighbor must be black 
P7. In an induced paw, the two odd-degree vertices must have different 

colors 

P8. In an induced C4, adjacent vertices must have different colors 
P9. If the neighborhood of any uncolored neighbor of a single (black) ver- 
tex s is contained in the neighborhood of s then the uncolored neighbor 
v of s minimizing weight(sv) must be black. If there are several options 
for vertex v, choose any one of them. We require rules PI and P8 to 
be applied before P9. 

Proof. The rules PI, P7, P8 follows from gj. while rules P3, P4, P5, P6 
follows from [SJ. 

The rule P2 follows from the coloring definition since each black ver- 
tex must be paired in order for the coloring to be valid. Finally, for P9, let 
s be a single vertex. Suppose the neighborhood of any uncolored neighbor 
of s lies within the neighborhood of s. Then the choice of the vertex to 
become the pair of s is independent of the choices for the remaining single 
vertices of the graph. Therefore, to obtain a minimum weighted dominat- 
ing induced matching of G, the neighbor v of s minimizing weight(sv) 
must be black, a 

Lemma 3. J3j/ If G contains a K4 then G has no DIM. 

Say that a coloring C is empty if all vertices are uncolored. 

Let C be a valid coloring and C an extension of it, obtained by the 
application of the propagation rules. If C = C' then C is called stable. On 
the other hand, if C 7^ C' then C' is not necessarily valid. Therefore, after 
applying iteratively the propagation rules, we reach an extension which 
is either invalid or stable. 
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In order to possibly extend a stable coloring C, we apply bifurcation 
rules. Any coloring directly obtained by these rules is not forced. Instead, 
in each of the these rules, there are two possibly conflicting alternatives 
leading to distinct extensions C[,C 2 of C. Each of C[ or C 2 may be 
independently valid or invalid. The next lemma describes the bifurcation 
rules. We remark that there exist simpler bifurcation rules. However, using 
the rules below we obtain a sufficient number of vertices that get forced 
colorings, through the propagation which follow the application of any 
bifurcation rule, so as to guarantee a decrease of the overall complexity 
of the algorithm. The complexity obtained relies heavily on this fact. 

In general, we adopt the following notation. If C is a stable coloring 
then S denotes the set of single vertices of it , U is the set of uncolored 
vertices and T = U\ U s< =sNu(s). 

Lemma 4. : BIFURCATION RULES 

Let C be a partial (valid) stable coloring of a graph G. At least one of the 
following alternatives can be applied to define extensions C[,C 2 of C . 

Bl. If C is an empty coloring: choose an arbitrary vertex v then C[ := 
C U BLACK{{v }) and C 2 := C U WHITE({v}) 

B2. If 3 edge vw s.t. v G N\j(s) and w G Nu(s'), for some s, s' G S, s ^ s' 
then C[:=CU BLACK({v}) and C 2 := C U WHITE({v}) 

B3. For some s G S, if3v£ Nu(s) s.t. 3w G N T (v): 

(a) // \Nu(s)\ / 3V d(w) / 3V \N T (v)\ > 2 then C[ := C U 
BLACK({v}) and C 2 := C U WHITE({v}). 

(b) = 3Ad(w) = 3AN T (v) = {w}, let Njj(s) = {v,v',v"}. 

i. // Nu(v') = Nu{v") = then C[ := C U BLACK({v}) and 

C' 2 :=C\J WHITE({v}) 

ii. // Nu(v') + 0, let w' £ N T (v'), with w' + w. If \N(w) U 
N(w')\ > 5 orww' <£ E(G) then C[ := C 'U BLACK ({v}) and 
C' 2 :=C\J WHITE({v}) 

iii. IfNjj{v') + 0, let w' £ N T (v'), with w' + w. If ww 1 G E(G) 
and z G N{w) n N(w') then C[ := C U BLACK({v"}), while 
if weight(sv) + weight(w' z) < weight(sv') + weight(wz) then 
C 2 :=CUBLACK({v}), otherwise C 2 := C ' \J BLACK \{v'}) 
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Each rule is applied after the previous rule, that is, if the condition of 
the previous case is not verified in the entire graph. Note that this applies 
to subitems of case B3. 

Proof If C is an empty coloring then the rule Bl is applied. 
If C is not an empty coloring and C is not a total coloring then S 7^ 0. 
Since C is not total and the graph is connected then there is at least one 
edge sv where v is uncolored. If s is white then v must be black P3else if 
s is a paired vertex then v must be white P4. Therefore s must be a single 
black vertex, hence S 7^ 0. Let s£S. Since C is valid then Njj{s) 7^ by 
V2 and since is stable \Nu(s)\ 7^ 1 by P6Therefore \Nu(s)\ > 2. Moreover 
rule P9can not be applied, therefore 3v G Nu(s) s.t. \Nu(v) \ N(s)\ > 0, 
let w G Nu(v) \ N(s). If 3s' G S,s ^ s' s.t. w G N v (s') then rule B2 is 
applied. 

Suppose that rule P2can not be applied. Then w G Nt{v){\Nt{v)\ > 1). 
Clearly, d(w) 7^ 1, otherwise, rule P2must be applied and v must get color 
black. 

In case \Nu(s)\ 7^ 3 or d(w) 7^ 3 or \Nt(v)\ > 2 we apply rule B3(a). 
Otherwise: \Nu(s)\ = 3,d(w) = 3, \N T (v)\ = 1. Note that in B3(b) when- 
ever we refer to v'w' it behaves symmetric to vw since otherwise v'w' were 
found in step B3(a) replacing vw. 

In the first subcase of B3(a) the case analyzed is whenever Njj{v') = 
Njj(v") = 0, while in the second and third the algorithm handle the cases 
when at least one of them has uncolored neighbors. 
Suppose w.l.o.g. Nu(v') 7^ where w' G Nj>{v'). It is easy to see that 
w 7^ w' since otherwise svwv' is a C4 and therefore w can't be uncolored 
by rule P8.Now there are three cases which lead to two possible outcomes 
from the algorithm: In case ww' G E(G) or \Nu(v) U Njj(w)\ > 5 then 
the result of the algorithm is given by the second subcase (ii), else it is 
given by the third subcase (iii). a 

4 The Algorithm 

The lemmas described in the last section lead to an exact algorithm for 
finding a minimum weight DIM of a graph G, if any, which we describe 
below. 

In the initial step of the algorithm, we find the set containing the K4S 
of G. If K4 ^ 0, by lemma [3j G does not have DIM's, and terminate the 
algorithm. Otherwise, define the set COLORINGS to contain through 
the process the candidates colorings to be examined and eventually ex- 
tended. This set should be implemented using a LIFO (Last In First 
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Out) data structure which achieves linear space complexity of the algo- 
rithm because the number of colorings in COLORINGS is at most n + 1, 
and each coloring needs 0(1) space. We give more detailed explanation in 
the next section. Next, include in COLORINGS an empty coloring. In the 
general step, we choose any coloring C from COLORINGS and remove it 
from this set. Then iteratively propagate the coloring by Lemma [2] into an 
extension C' of it, and validate the extension by Lemma [TJ The iterations 
are repeated until one of the following situations is reached: C' is invalid, 
C is a total valid coloring, or a partial stable (valid) coloring. In the first 
alternative, C is discarded and a new coloring from COLORINGS is 
chosen. If C total valid coloring, then sum the amount of valid 

DIMs related to this coloring, find its weight and if smaller than the least 
weight so far obtained, it becomes the current candidate for the mini- 
mum weight of a DIM if G. Finally, when C is stable we extended it by 
bifurcation rules: choose the first rule of Lemma [5] satisfying C , compute 
the extensions C' and C" , insert them in COLORINGS , select a new 
coloring from COLORINGS and repeat the process. 

The formulation below describes the details of the method. The propa- 
gation and validation of a coloring C is done by the procedure PROPAGATE— 
V ALI DAT E(C, RESULT). At the end, the returned coloring corre- 
sponds to the extension C' of C, after iteratively applying propagation. 
The variable RESULT indicates the outcome of the validation analysis. If 
C' is invalid then RESULT is returned as 'invalid'; if C is a valid total 
coloring then it contains 'total', and otherwise RESULT equals 'partial'. 
Finally, BIFU RCATE(C, C[,C' 2 ) computes the extensions C[ and C' 2 of 
C. 
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Algorithm Minimum Weighted DIM / Counting DIM 



1. 


Find the subset KA 




if KA ^ then terminate the algorithm: G has no DIM 






SOLUTION := NODIM 


2. 


COLORINGS := {C}, where C is an empty coloring 


o 


while COLORINGS ^ do 




a. 


choose C e COLORINGS and remove it from COLORINGS 




b. 


PROPAGATE - VALID ATE(C, RESULT) 




c. 


if RESULT = 'total' and weightiC) < SOLUTION then 






SOLUTION := weight(C) 






else if RESULT = 'partial' then 






Set C( and C£ according to BIFURCATION RULES on C Bl 






COLORINGS := COLORINGS U {C{, C^} 






end if 


4. 


Output SOLUTION 



procedure PROPAGATE - VALID AT E(C, RESULT) 
Comment Phase 1: Propagation 

1. a := c 

2. repeat 

c := a 

C := extension of C obtained by the PROPAGATION RULES 
until C = C 
Comment Phase 2: Validation 

3. Using the VALIDATION RULES (lemma Q} do as follows: 
if C is an invalid coloring then return (C, 'invalid') 
else if C is a partial coloring then return (C, 'partial') 
else return (C, 'total') 



5 Correctness and Complexity 

It is easy to see that our algorithm uses the branch & reduce paradigm 
since propagation rules can be mapped to reduction rules since are used 
to simplify the problem instance or halt the algorithm and the bifurcation 
rules can be mapped to branching rules since are used to solve the problem 
instance by recursively solving smaller instances of the problem. 
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Theorem 1. The algorithm described in the previous section correctly 
computes the minimum weight of a dominating induced matching of a 
graph G. 

Proof: The correctness of the algorithm follows from the fact that 
the algorithm considers all the cases that need to be considered, this is, 
any coloring that represents a DIM must be explored. Lemmas [2] and 2] 
ensures that the simplifications of the instances are valid, invalid colorings 
are discarded, some valid colorings can not be explored only if other valid 
coloring representing a better DIM (with less weight) is explored. 

For proving the worst-case running time upperbound for the algorithm 
we will use the following useful definition and theorem. 

Definition 2. flSf Let b a branching rule and n the size of the instance. 
Suppose rule b branches the current instance into r > 2 instances of size 
at most n — t\, n — t2, ■ ■ ■ , n — t r , for all instances of size n > max{ti : 
i = 1,2, .. . , r}. Then we call b = (ti,tz, . . . ,t r ) the branching vector of 
branching rule b. 

The branching vector b = (t\, t2, ■ ■ ■ , t r ) implies the linear recurrence 
T{n) < T{n - t x ) + T(n - t 2 ) + . . . , T(n - t r ). 

Theorem 2. jll3f Let b be a branching rule with branching vector (ti, t2, ■ ■ ■ ,t r ). 
Then the running time of the branching algorithm using only branching 
rule b is O*(o n ) ; where a is the unique positive real root of 

rj.n rj,n—t\ y.n—t2 ^,n—t r _ q 

The unique positive real root a is the branching factor of the branching 
vector b. 

We denote the branching factor of (ti, t2, ■ ■ ■ ,t r ) by r(ti,t2> • • • > tr)' 
Therefore for analyzing the running time of a branching algorithm we 
can compute the factor a, for every branch rule bi, and an upper bound 
of the running time of the branching algorithm is obtained by taking 
a = maxiOti and the result is an upper bound for the running time of 
O*(o n ). 

The upper bound comes from counting the leaves of the search tree 
given by the algorithm, using the fact that each leave can be executed 
in polynomial time. The complexity of the algorithm without hiding the 
polynomial depends on the upperbound time for the execution of each 
branch in the search tree. 

Further notes on this topic can be found in |13| 
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Theorem 3. The algorithm above described requires 0*(1.1939 n ) time 
and 0(n + m) space for completion. 

Proof: 

Using the definition [2] and the theorem [2] the calculation of the up- 
per bound time is reduced to calculation of the branching vector for each 
branching rule (i.e. bifurcation rules in our algorithm) and obtain the 
associated branching factor for each case. Then the bound is given by 
the maximum branching factor. Note that to use this we must observe 
that the reduction rules (i.e. propagation rules in our algorithm) can be 
computed in polynomial time and leads to at most one valid extension 
of the considered coloring. So, the propagation rules do not affect the 
exponential factor of the algorithm. Moreover, each branch of the algo- 
rithm has cost 0{n + m) in time and space. This is easy to note since 
from the empty coloring up to any total coloring each vertex v is painted 
once and the cost in time incurred for painting each vertex is given by 
the updating of the color of the vertex and updating this information for 
the neighborhood, hence |iV(t?)| times a constant operation for updating 
a counter with amount of black/white/uncolored neighbors. Therefore, 
the total cost for each branch is 0{n + m). 

Let's analyse each bifurcation rule to obtain the maximum branching 
factor. 

1. If C is an empty coloring: choose an arbitrary vertex v then C[ := 
C U BLACK{{v}) and C' 2 := C U WHITE{{v}): It is easy to see 
that this rule is executed once, after that, the coloring is never empty 
again. Since this rule bifurcation opens two branchs then we can up- 
per bound the time of the algorithm by 2 times the complexity of 
the algorithm executed in an instance of size n — 1. Therefore the 
asymptotic behavior of the algorithm is not affected. 

2. If 3 edge vw s.t. v G Njj{s) and v' € Nu(s'), for some s,s' G S,s ^ s' 
then C[ := C U BLACK({v}) and C 2 := C U WHITE({v}). 

Here we extend the original coloring C to C[ and C' 2 by coloring the 
vertex v with black and white respectively. Recall that exists an edge 
vw such that v £ Nu(s), w £ Nu(s'). If v is black then Nu(s) \ v are 
white, while v' is white. On the other hand, if v is white then w is 
black and Njj(s')\w are white. Therefore the size of uncolored vertices 
is reduced for each branch (i.e. for each new coloring). The associated 
branching vector is (1+ |i%(s)|, 1+ \Nu(s')\). By rule P2 \N v {s)\ > 2 
and \Njj(s')\ > 2. The following observation turns out to be useful: 
If \Njj(si)\ = 2 then Njj(si) can be totally painted wether v is black 
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or white. The case with N u (s') is symmetric. Therefore the branching 
vector with biggest branching factor is (3,5) (r(3, 5) « 1.1939). 
3. For some s G S, if 3w G i%(s) s.t. 3w G N T (v): 

Note that if /Siu G Nt(v) for any i> G Njj(s) then either the propagat- 
ing rule P9or P5can be applied to get an extension of the coloring. 

(a) If \Nu(s)\ / 3V d(w) / 3 V \N T (v)\ > 2 then C[ := C U 
BLACK({v}) and C 2 := C U JT.E({v}). 

Since w is uncolored then w is not a pendant vertex, d(w) > 1. 
Since U; is uncolored then it has nor white nor paired black neigh- 
bor. Moreover, if w has a single black neighbor then this is the 
case analyzed above. Therefore w has uncolored neighbors and let 
x be one of them. 

(a.l) \N T (v)\ > 2: Let v' G N T (v). In C[ {v,x} will be black 
while {vi,v',w} will be white. In C 2 {v} will be white while 
{v',w} will be black. This lead to the branching vector (3,5). 

(a.2) d(w) ^ 3. Hd(w) = 2 then in C[ the vertices Nu(s)U{w, x} 
will be colored and in C' 2 the vertices {v, x} will be black while 
{w} will be white. Therefore the branching vector will be at 
least (3,5). 

Else if d(w) > 3 then in C[ the vertices Nu(s) U Njj[w] will be 
colored and in C' 2 the vertices {v, w} will be colored. In case 
\Nu{s)\ = 2 then v\ will be colored too. Therefore the branch- 
ing vector (2,7) (t(2,7) = 1.1908). 

(a.3) \N v (s)\ = 2: Let N v {s) = and N(w) = {v,x,x'}. 

In C[ after applying propagation rules the vertices {v, x, x'} 
will be black while {v±,w} will be white. In C 2 after applying 
propagation rules the vertices {v±,w} will be black while {v} 
will be white. The result is the branching vector (3,5). 

(a.4) \Nu(s)\ > 3: Let {vi,v 2 ,v 3 } G Nu(s) and N(w) = {v,x,x'}. 
In C[ after applying propagation rules the vertices {v, x, x'} 
will be black while {v±, v 2 , vs, w} will be white. In C' 2 after ap- 
plying propagation rules the vertices {w} will be black while 
{v} will be white. The result is the branching vector (2,7) 
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(b) // \N v (s)\ = 3 A d(w) = 3 where Nu(w) = {v, x, x'}, N v (s) = 
{v, v' , v"} Note that {x, x'}n{v, v 1 , v"} = since otherwise at least 
one of them must be colored by rule P8. 
(b.l) IfNu(v') = Nu(v") = then 

C[:=CU BLACK{{v}) and C' 2 :=C\J WHITE{{v}) : 
Suppose w.l.o.g. weight(sv') < weight(sv"), then 
In C[ after applying propagation rules the vertices {v,w',w"} 
will be black while {v ' , v" , w} will be white. In C 2 after apply- 
ing propagation rules the vertices {v',w} will be black while 
{v,v"} will be white. The result is the branching vector (4,6) 
(r(4,6) = 1.1510). 

(b.2) IfNu(v') + 0, let w' G N T {v'), with w' + w: 
If \N T [w] U N T [w']\ > 5 then 

C[ := C U BLACK({v}) and C 2 := C U WHITE({v}) 

Note that if d(w') ^ 3 then v'w' satisfies the properties of an 
already analized case, hence C[ := C U BLACK({v'}) and 
C 2 := CUWHITE({v'}). 

Since d(w) = d(w') = 3 and |./Vt[w;] U A^[it; / ]| > 5, then 3x,y 
s.t. x G N T (w),x N T (w') and y G N T (w'),y £ N T (w). In 
C[ after applying propagation rules the vertices {v, x, x' , w'} 
will be black while {v',v",w} will be white. If x' = w' then y 
must be black by rule P6.In C' 2 the vertex {w} will be black 
while the vertex {v } will be white. The result is the branching 
vector (2,7) 

(b.3) If Nu{v') / 0, let w' G N T {v'), w' + w 

If \N T [w] U N T [w']\ < 3 and z£ N(w) n N(w') then 
C[ :=CUBLACK({v"}), 

if weight(sv) + weight(w' z) < weight(sv') + weight(wz) then 
C' 2 :=CUBLACK({v}) 
otherwise C' 2 := C U BLACK{{v'}) 

Since d(w) = d(w') = 3 then ww' G -E(G) and 3z G Nt(v) n 
Nt(w), otherwise the case is one of the above. 
In both colorings, C[ and C2 the vertices {v, v' , i/', t/;', z} will 
be colored. The branching vector is (6,6). (r(6, 6) = 1.1225). 

The worst branching factor is r(3, 5) ~ 1.1939. In consequence, the 
time complexity of this algorithm is O * (1.1939 ra ). 
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To achieve linear space complexity, we use a stack to store the color- 
ing sequence of the current branch. The only additional space is needed 
for COLORINGS and extra information to restore the initial condition 
for each coloring. For each coloring c € COLORINGS extended from a 
bifurcation rule we store the number of colored vertices before the bifur- 
cation, the vertex colored during bifurcation and its color. These elements 
are sufficient to restore the initial condition, a 

The analysis can be extended for the case of non-connected graphs. It 
is easy to obtain the same upper bound after separating the cases where 
each connected component of four or less vertices is solved in constant 
time. 

6 Counting the number of DIM's 

The previous algorithm can be easily adapted to count the number of 
DIM's. The number of DIM's is the number of total valid colorings. Given 
a coloring C we define TVC{C) the number of total valid colorings that 
can be extended from C. If we apply any propagation rule to coloring C 
we obtain a coloring C' . Clearly TVC{C) = TVC(C), except for rule 
P9. In the later case TVC(C) = TVC{C) ■ \N v {s)\ where s is the single 
vertex chosen to apply the rule. 

If we apply any bifurcation rule to coloring C we obtain two extended 
colorings C[ and C' 2 . Clearly TVC{C) = TVC(C[) + TVC{C 2 ), except 
for rule B3(b)iii. In the later case TVC{C) = TVC(C[) + 2 • TVC{C 2 ). 

Using the above facts it is trivial to modify the algorithm to solve the 
counting problem. 

7 Conclusions 

We have developed a new exact exponential algorithm for an extensively 
studied problem. Moreover the developed algorithm is practical since 
there are no big constants or polynomials hidden in the upper-bound 
and it is straightforward to implement it. 



Problem 


Previous results 


New results 


Weighted DIM 


0*(1.4423 m ) [4 21 231, 0(1.4423" • m) [17 \ 


0(1.1939" -m) 


Counting DIM 


0*(1.3247 m ) [4121110] . 0(1.4423" • m) p2] 


0(1.1939" -m) 
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